A new approach to reconstruction methods in f(R) gravity 
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We present a new approach of the reconstruction method based on the use of the cosmic param- 
eters instead of a time law for the scale factor. This allows the derivation and analysis of a set of 
new non-trivial cosmological solutions for /(i?)-gravity. A number of simple examples are given. 
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INTRODUCTION 

The ACDM (or Concordance) Model [l| which is based 
on the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric is one of the greatest successes of General Rel- 
ativity. It reproduces beautifully all the main observa- 
tional results e.g., the dimming of type la Supernovae 
.2], Cosmic Microwave Background Radiation (CMBR) 
anisotropies [3(, Large Scale Structure formation [4j, 
baryon oscillations [5j and weak lensing [6|). Unfortu- 
nately, this model is also affected by significant fine- 
tuning problems related to the vacuum energy scale and 
therefore it is important to investigate other viable the- 
oretical schemes compatible with cosmological observa- 
tions. 

A now widely studied alternative to the ACDM model 
is based on gravitational actions which are non-linear in 
the Ricci curvature R and/or contain terms involving 
combinations of derivatives of R: the f(R) theories of 
gravity [9l4l2|. Such models first became popular in the 
1980's because it was shown that they are derived from 
fundamental physical theories (Like M-theory) and natu- 
rally admit a phase of accelerated expansion which could 
be associated with an early universe inflationary phase 
[8J]. The fact that the phenomenology of Dark Energy 
requires the presence of a similar phase (although only a 
late time - low energy one) has recently revived interest 
in these theories. In particular, the idea that Dark En- 
ergy may have a geometrical origin, i.e., that there is a 
connection between Dark Energy and a non-standard be- 
havior of gravitation on cosmological scales has resulted 
in it becoming a very active area of research over the past 
few years (see for example [13l4l8|). 

Unfortunately efforts to obtain an understanding of the 
physics of these theories are hampered by the complexity 
of the fourth-order field equations, making it difficult to 
obtain both exact and numerical solutions, which can be 
compared with observations. Recently, however, progress 
has been made in resolving this issue using a number of 



useful techniques. One such method, based on the theory 
of dynamical systems [19j | has proven to be very successful 
in providing a simple way of obtaining exact solutions 
and a (qualitative) description of the global dynamics of 
these models |20| . 

Another interes ting technique is the so called recon- 
struction method 21|. Here one assumes that the ex- 



pansion history of the universe is known exactly and one 
inverts the field equations to deduce what class of f(R) 
theories give rise to a given FLRW model. The exis- 
tence of such solutions is particularly relevant because in 
FLRW backgrounds they typically represent asymptotic 
or intermediate states in the full phase-space of all pos- 
sible cosmological evolutions. In |22| it was found that 
expansion histories based on a power-law solution for the 
scale factor and a perfect fluid only exist for R n gravity, 
while in [23| it was shown that in order to reconstruct 
a A CDM expansion history in f(R) gravity, it is neces- 
sary to add additional degrees of freedom to the matter 
sector. 

A fundamental limitation of the reconstruction method 
is that, because of the mathematical steps involved, only 
very simple cosmic histories (i.e., single power law behav- 
iors) can be successfully connected to a /(-R)-theory in an 
exact way [24|. On the other hand /(i?)-theories of grav- 
ity, being of order four in general admit a much richer 
set of solutions, which have at most four different modes. 
From this point of view it is clear how reconstructing 
with a single power law behavior corresponds to a dra- 
matic constraint on the form of the function /. As a 
consequence one obtains by reconstruction, Lagrangians 
which arc often of little use due to their complexity. 

In this paper we propose a modification of the clas- 
sic reconstruction scheme able to generate more general 
cosmic histories and able to give much simpler results 
for the form of f(R). This will be done describing these 
expansion histories by means of the cosmic parameters 
rather than a function of time. 



FIELD EQUATIONS FOR f(R) FLRW MODELS 

Let us consider the typical action for /(i?)-gravity in 
the c = 8ttG = 1 units: 



A = / <PxJ~g 



f(R) + C, 



(1) 



where R is the Ricci scalar, / is general differentiable (at 
least C 2 ) function of the Ricci scalar and C m corresponds 
to the matter Lagrangian. 

In a FLRW universe, the field equations take the form 



3H + 3H 2 


1 

-~2f 


p + 3p + f- f'R + 3Hf"R 


+3f"R 2 + 3f"R 


(2) 


™' = 7 


p+ RfLJ.-3Hf»k 


(3) 



i.e., the Raychaudhuri and Friedmann equations. Here H 
is the Hubble parameter, which defines the scale factor 
a(t) via the standard relation H = a/a, the Ricci scalar 



R = 6H+12H 2 



(4) 



and /',/" and /'" abbreviates d n f/(dR) n for n = 1..3 
respectively. The energy conservation equation for stan- 
dard matter 



p=-3H(p + p) 



(5) 



closes the system. 

It is interesting to note that the Raychaudhury equa- 
tion can be obtained by adding the Friedmann equation 
to its time derivative and using the Energy conservation 
equation and the definition of the Ricci scalar. Hence, 
any solution of the Friedmann equation automatically 
solves the Raychaudhuri equation. Thus, in the recon- 
struction process we only need to solve the Friedmann 
equation, as that guarantees a solution to the other equa- 
tions. 



RECONSTRUCTION STRATEGIES 

In this section we introduce, after reviewing briefly the 
standard reconstruction method, a new way of obtain- 
ing the form of the function f(R) starting from condi- 
tions on the dynamics of the cosmic parameters. We also 
give some specific examples with each of these strategies, 
which shows that this can be a powerful tool for gener- 
ating exact solutions in fourth order gravity. 

The classical reconstruction method 



In the classical reconstruction method 



21(, the func- 
tion f(R) is derived from a given solution a = a(t). An 



explicit reconstruction is possible if and only if the func- 
tion R(t), as obtained by substituting the solution in 
equation (|4|, is invertible analytically. That is, if we 
have an explicit function g such that t = g(R). When 
this is possible, we can write all the functions, (namely 
a(t), H(t), R(t)) in terms of the variable R and the Fried- 
mann equation becomes 

3H[g(R)]R[g(R)} f" + f (3H[g(R)} 2 - R 



I 
2 



Pog{R) 



-3(1+™) 



(0) 



where we have assumed the matter to be a perfect fluid 
with equation of state p = wp. The above equation is 
now a second order differential equation for the function 
f(R), the solution of which gives the class of theories of 
gravity for which the given function a — a(t) is an exact 
solution. 

As an example let us consider the power law solution 
a(t) = aot m [22J], From Q we see that the Ricci scalar 
is invertible and is given by 



R = 6m(2m- l)t 



- 2 -- a t~ 2 



(7) 



Solving the Friedmann equation we obtain the following 
general solution [22| 



f(R) = A mw [ — 



ym(l+1f) 



Ci-RI" 



1% 



,C 2 RT 



(8) 



where A mw and j3 m are constants depending on the values 
of m and w and po, and C\ t % are arbitrary constants of 
integration. This is a lucky coincidence. Considering 
more general solutions leads at best to a combination of 
hypergeometric functions (see e.g. 2l|). 



Reconstruction from the condition a — h(a) 

Let us consider now consider the relation d = h(a) 
instead of a = a(t). This equation relates the Hubble 
parameter to the scale factor. We can easily calculate 
the Ricci scalar in this case as 

'l(/l 2 ),q h 2 ' 

2 



R(a) = 6 



(9) 



Like the previous case, an explicit reconstruction is pos- 
sible if the function R(a), is invertible analytically, that 
is, if we have an explicit function g such that a = g(R). 
However, this inversion is now possible in cases which 
were forbidden before. 

The Friedmann equation can now be written as 

o %(fl)] p r (m]h , fm]f/l , fl \. Ag{R)] 2 Rf 

3 ^(kr R ' a[9{ R)] %( R)]f +/ r^r + ~ 



2 



P0S(i*)- 3(1+,o) . 



(10) 



The solutions to the above equation gives the class of 
theories of gravity for which the condition a = h(a) is 
satisfied. This last relation can be then integrated to 
obtain the corresponding scale factor : 



t = 



a da 
o h (a) 



(11) 



As an example let us investigate for which class of model 
a dust-like matter (w = 0) behaves in the following way 



• m i — r^ 

a = — = va- Aa z 
VA 



(12) 



As one can easily see, this condition corresponds to a 
cyclic universe with the scale factor 



a(t) = ^sm 2 (nt) 

VA 



(13) 



Using equation (0) , we get an analytically invertible Ricci 
scalar 



R(a) = 



12ft 2 (3 - 4Aa) 
Aa 



(14) 



Substituting everything in equation (|10[) , we find the par- 
ticular solution to be 

f{R) = a 1 R + a 2 R 2 + a 3 R 3 + a 4 , (15) 

where a n {n — 1..4) are constants depending on A, Q and 
Po- 



Solving the above equation we obtain the f(R) theory 
which admits a solution compatible with the condition 
above. To solve for the scale factor as a function of time 
we express this condition as a differential equation for 
a(t): 



, a a \ , ,a. 

6 - =h(-) 

a a z I a 



(19) 



Let us illustrate this with a simple example by assuming 
a vacuum (p = 0) universe and let the condition be 



H = m 



(20) 



where m is a constant. In this case we get a Ricci scalar 
which is invertible and the solution of the Friedmann 
equation gives the following theory of gravity: 



f(R)=a 1 R + a 2 R 2 +a 3 



(21) 



where the constants a n (n — 1..3) depends on m. Now 
solving for the scale factor we obtain 



a(t) = a exp (y (t 2 - C x i)) , 



(22) 



which, for m > 0, represents a universe bouncing in the 
past. 

Note that this solution, unlike the previous ones con- 
tains two integration constants. This is the consequence 
of the fact that H = h(H) is a second order differential 
equation in a [25[. 



Reconstruction from the condition H = h(H). 

Following the line of reasoning above, one can also 
imagine describing the scale factor as a differential equa- 
tion for the parameter H e.g., H = h(H). In this case, 
we can immediately find the Ricci Scalar in terms of the 
Hubble parameter: 



R{H) = 6h(H) + 12H 2 



(16) 



and an explicit reconstruction is possible if the above 
relation is analytically invertible, that is, we have an ex- 
plicit function g such that H — g(R). The scale factor in 
terms of the Ricci scalar can then be found solving the 
integral 



a(R) — exp 



g(R)dg(R) 

h(g(Rj) 



(17) 



Substituting for the above quantities, the Friedmann 
equation becomes 



3g(R)R, H \g(R)]h\g(R)]f" + f 3g(R) 



Rf 



f 



1^ 



a (£)3(l+«0 



(18) 



Reconstruction from the condition q — h(q) 

It is now clear how we can further generalize the above 
strategy. One can now give a condition on the dynamics 
of the declaration parameter q = — da/ a 2 . 

Supposing q = h{q), we can integrate to find the Hub- 
ble parameter in terms of q as 



1 + 9 



dq 



(23) 



H{q) J h{q) 

The Ricci scalar in terms of q is given by 

R(q) = 6H(q) 2 (l~q). (24) 

For explicit reconstruction, the above equation should be 
analytically invertible as in other cases, so that q = g(R). 
We can then solve for the scale factor in terms of the Ricci 
scalar as 



*-(/m*< b ' 



h(g(R)) 

Substituting into the Friedmann equation, we obtain 



(25) 



3H[g(R)}R, q [g(R)}h[g(R)}f 



2 Rf] , / 



-f{ 3H(g(R)) 2 -j- 2 



Po 



a (R)m+™) 



,(26) 



the solution of which gives the required theory. 
As an illustration let us assume 

q = to(1 + q)y/q , 



(27) 



where m is a constant. Reconstructing the theory for a 
universe filled with dust-like matter, we get the following 
function f(R): 



f(R) = aRy/R - 3m 2 



(28) 



where a is a constant depending on p$. Using equation 
([24)) . we can then solve for the scale factor as a function 
of time, and is given by 

a{t) = d sin ( —mt J + C 2 cos ( —mt J , (29) 

which has oscillatory behaviour. Note that the solution 
above has, as expected, three constants due to the fact 
that q = h(q) is a third order differential equation in a. 



Reconstruction from the dynamics of higher order 
parameters 

Since the f(R) theories are of order four it is natural 
to give a relation for the scale factor involving a fourth 
order equation for a. This would involve naturally the 
jolt parameter j = da 2 /a? \J\. It turns out that with 
this parameter makes it difficult to achieve the results of 
the previous section. 

For this reason, in analogy with what is usually done 
with the scale factor, we propose the following expansion 
of the Hubble parameter: 



H = H + H(t- t ) + -H(t - to) 2 + 



(30) 



which allows us to define two new cosmic parameters: 



*-■§• 



3 = 



HH 



(31) 



igffi, so that 



It is clear that Q = —(q + 1) and J 
they are related to the standard parameters j and q 
We can easily check from the above definitions that 



1 Q 



= 3-2. 



(32) 



HQ 2 

Let us now suppose that we would like to reconstruct 
from the given condition 

1 



H 



3 = h{J) 



(33) 



We then integrate the above conditions to find 



1 



Q{3) 



2-3 



d3 ; H(3) = exp 



0(3) 
h{3) 



dJ 



(34) 



a(J) — exp 



dJ 

W) 



(35) 



The Ricci scalar can now be written in terms of these 
quantities as 



i?(J) = 6#(J) 2 (2 + Q(J)). 



(36) 



If the above relation is explicitly invertible, that is we 
have J — g(R), then we can again substitute everything 
into the Friedmann equation to give 



/' ^ 3H(g(R)) 



3H[g(R)}Rj[g(R)}H[g(R)}h[g(R)]f" 
2 Rf'\ , / Po 



a[ 5 (i?)] 3 ( 1 +-) 



0(37) 



Solving this equation gives the required theory. 

As a specific example, let us assume the condition to 
be 



1 
H' 



-3 = 3-2 



(38) 



Also for simplicity let us consider a universe filled with 
dust-like matter. Solving the Friedmann equation we get 
the following particular solution: 



f(R) = aii? + a 2 R 2 + a 3 R 3 



«4 



(39) 



where a n (n = 1..4) are constants depending on po- Now 
solving for a(t), we obtain the following solution in inte- 
gral form: 



a(t) 



d.v 



yV + 2x 3 + CI 



(40) 



It is clear that this solution, although given implicitly, 
contains four different modes. 



DISCUSSION AND CONCLUSION 

In this paper we have presented a new approach to re- 
construction methods for /(i?)-gravity. Using the stan- 
dard and some newly defined cosmic parameters, we have 
been able to connect some non-trivial cosmic histories to 
some relatively simple f(R) functions, solutions which 
would have been harder to obtain through a direct in- 
tegration of the field equations. Given the simplicity of 
the functions f(R) obtained, one can proceed to examine 
other features of these solutions such as the evolution of 
cosmological perturbations in these models. 

The few examples presented here are an indication of 
the richness of the general behavior of /(i?)-cosmology, 
which so far could be grasped only via the dynamical sys- 
tem approach. We are confident that this new method 
will help uncover more details of the cosmology of these 
theories and lead to a deeper understanding of their fea- 
tures and structure. 
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